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Abstract 

I We present a fully microscopic three-cluster nuclear model for light nuclei on the 

' basis of a J-Matrix approach. We apply the Modified J-Matrix method on ^He and ^Be 

. for both scattering and reaction problems, analyse the Modified J-Matrix calculation, 

and compare the results to experimental data. 

o 

1 Introduction 

^ , The J-Matrix method (JM) has proven very successful in microscopic nuclear calculations, 

pH I particularly the Modified J-Matrix approach (MJM)[lj^, also called the Algebraic Model 
in some nuclear physics literature. Both collective [3] and cluster descriptions |11E|, of light 
(p-shell) nuclei have been studied with the MJM, using an oscillator basis. 

The Harmonic Oscillator basis has always been very popular in nuclear physics. For light 
nuclei, spherical oscillator states have often been used as a first approximation to the single- 
particle orbital wave functions in the popular nuclear shell-model. The nuclear many-body 
basis is then built up as a set of Slater determinants of single-particle oscillator orbital states 
to take the Pauli principle into account. 

Many nuclear two-body potentials feature a superposition of Gaussian components. Ma- 
trix elements of two-body operators for Slater determinants reduce to a simple sum of two- 
body matrix elements, involving the single-particle oscillator states. A Gaussian form of the 
operator then easily leads to an analytical form for the matrix elements. This immediately 
shows the computational advantage of using an oscillator state (or even a superposition of 
oscillator states) for the single-particle wave functions in a fully microscopic many-body 
nuclear model. 
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One of the features of the oscillator basis is the Jacobi (tridiagonal) form of the matrix 
of the kinetic energy operator. This makes it a proper candidate for considering the JM 
approach to solve the Schrodinger equation expressed in matrix form. 

The nuclear many-body system is very complex though, and requires huge superpositions 
of shell-model states to reproduce spectral properties over an important energy range. Also, 
the nuclear system exhibits several modes when the system is energetically excited, and 
possibly fragmented. There is an interplay between collective modes, such as monopole 
and quadrupole excitations, and cluster effects that are particularly pronounced when the 
nucleus disintegrates. One therefore often introduces very specific antisymmetrized forms for 
the wave function with a specific configuration of single-particle orbitals, that feature some 
specific collective behavior that one wants to study. In this way the Hilbert space is limited 
to a single (or a few coupled) nuclear model state (s) in which the collective coordinates are 
the only remaining dynamical coordinates. The dimensions of the Schrodinger equation are 
then strongly reduced. 

Although the nuclear two-body interaction has a short range, the effective interaction as 
a function of the collective coordinates often displays a long range. If the collective behavior 
is described as a superposition of oscillator states, this is usually also reflected in the matrix 
elements of the nuclear potential, which display a slow decrease for increasing oscillator 
excitation. This clearly limits the applicability of the standard JM method, as too large 
energy matrices have to be considered in the internal, non-asymptotic, region. Indeed, the 
main computational cost for nuclear JM calculations often lies in the construction of the 
matrix elements of highly excited states. To account for this limitation the MJM approach 
was developed ^113, introducing a semi-classical approximation for the matrix elements in 
the highly excited internal region. This leads to a modification of the standard three-term 
recursion relation for the far interaction and asymptotic region by including (semiclassical) 
potential contributions. A drastic reduction of the matching position for boundary condition 
is obtained, resulting in a remaining matrix equation of manageable dimensions. 

The Coulomb contribution to the nuclear potential, which is known to have a long range, 
can be handled in the same way. 

In the following sections we will elaborate on the cluster description of nuclei as an 
important collective description for the lightest nuclei, in which the internuclear distances 
then become the dynamical (collective) coordinates. We will link this model to the MJM 
approach to solve the Schrodinger equation, and present applications on 6-particle nuclear 
systems where three-cluster effects are important. We will compare the theoretical results 
to the experimental ones to test the validity of the approach. 

2 JM Cluster models for light nuclei 

Two- and three-cluster configurations carry an important part of the low-energy physics 
in light nuclei, and directly relate to scattering and reaction experiments in which smaller 



2 



N 




N 




Figure 1: Two configurations of Jacobi coordinates for the three-cluster system a + N + N . 



fragments are used to study compound properties. 

In this section we discuss the three-cluster description for light nuclei. Where appropriate 
we briefly present some two-cluster properties, as the three-cluster approach is essentially a 
generalization hereof. 

The many-particle wave functions for a three-cluster system of A nucleons {A = Ai + 
A2 + A^) can be written, using the anti-symmetrization operator A, as follows 

^ (qi, .., q^_i) = A (A) ^2 {A2) ^3 {A3) [R)] (1) 

where the centre of mass of the A-nucleon system has been eliminated by the use of Jacobi 
coordinates so that only internal dynamics are described. The cluster wave functions 

^^ (A,) 

vi/,(A,) = M/,(q«,..,qg_i) = 1,2,3) (2) 

represent the internal structure of the i-th cluster, centered around its centre of mass Rj. To 
limit the computational complexity of the problem, these cluster functions are fixed and they 
are Slater determinants of harmonic oscillator (Os)-states, corresponding to the groundstate 
shell-model configuration of the cluster {Ai < 4 for all i). The \E'ij (R) wave function 

vj/^ (R) = qf ^) = ^i? (qi, q2) (3) 

represents the relative motion of the three clusters with respect to one another, and qi and 
q2 represent Jacobi coordinates. In figure ^ we indicate an enumeration of possible Jacobi 
coordinates and their relation to the component clusters. 

The state (jH)) is not limited to any particular type of orbital; on the contrary we will 
use a complete basis of harmonic oscillator states for the relative motion degrees of freedom. 
Thus the full A-particle state cannot be expressed as a single Slater determinant of single 
particle orbitals. 
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An important approximation, known as the "Folding" model, is obtained by breaking the 
Pauli principle between the individual clusters, but retaining a proper quantum-mechanical 
description of the clusters, which is described by the wave function 

(qi, ., q^_i) = ^1 (Ai) ^2 {A2) ^3 (^3) ^iJ {R) (4) 

Because each cluster wave function is antisymmetric (they are Slater determinants) one is 
indeed neglecting the inter-cluster anti-symmetrization only. The Folding Model has the 
advantage of preserving the identities of the clusters and, if the intra-cluster structure is 
kept "frozen" , it reduces the many-particle problem to that of the relative motion of the 
clusters. 

For a two-cluster description the formulae simplify to: 

m (qi, .., q^-i) = A [*i (Ai) ^2 {A2) {R)] (5) 

with 

vl/H(i?) = v&^(qi''^) =^i?(qo) (6) 

The folding approximation will be the natural choice for calculating the asymptotic be- 
havior of the cluster-system, i.e. the disintegration of the system in the two or three non- 
interacting individual clusters. This amounts to the situation that all clusters are a sufficient 
distance apart and inter-cluster antisymmetrization has a negligible effect. 

Because the cluster states are fixed and built up of (Os)-orbitals, the problem of labeling 
the basis states with quantum numbers relates to the inter-cluster wave function only. This 
holds true whether one uses the full anti-symmetrization or the folding approximation. In a 
two-cluster case, the set of quantum numbers describing inter-cluster motion is unambigu- 
ously defined, and is obtained from the reduction of the symmetry group f/(3) D 0(3) of the 
one-dimensional oscillator. This reduction provides the quantum numbers n for the radial 
excitation, and L, M for the angular momentum of the two-cluster system. 

In a three-cluster case, several schemes can be used to classify the inter-cluster wave 
function in the oscillator representation. In [HI Hj three distinct but equivalent schemes 
were considered. One of these used the quantum numbers provided by the Hyperspherical 
Harmonics (HH) method (see for instance [HI El Ell)- This is the classification that we will 
adopt. Even within this particular scheme there are several ways to classify the basis states. 
We shall restrict ourselves to the so-called Zernike-Brinkman basis This corresponds to 
the following reduction of the unitary group f/(6), the symmetry group of the three-particle 
oscillator Hamiltonian, 

f/(6) dO(6) dO(3)®0(3) dO(3) (7) 

This reduction provides the quantum numbers K, the hypermomentum, n, the hyperradial 
excitation, /i, the angular momentum connected with the first Jacobi vector, /2, the angu- 
lar momentum connected with the second Jacobi vector, and L and M the total angular 
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momentum obtained from coupling the partial angular momenta /i, l2- Collectively these 
quantum numbers will be denoted by v, i.e. u — {n, K, {lil2)LM} in the remainder of the 
text. 

There are a number of relations and constraints on these quantum numbers: 

• the total angular momentum is the vector sum of the partial angular momenta hand 
I2, i.e. L = I1+ I2 or |/i -kl < L <li + I2. 

• by fixing the values of /i and I2, we impose restrictions on the hypermomentum K = 
h + h, /i + /2 + 2, /i + /2 + 4, . . . This condition implies that for certain values of 
hypermomentum K the sum of partial angular momenta l\ + I2 cannot exceed K. 

• the partial angular momenta li and I2 define the parity of the three-cluster state by 
the relation 7r= 

• for the "normal" parity states tt = (—1)'^ the minimal value of hypermomentum is 
i^'min = whereas i^Tmin = L + 1 for the so-called "abnormal" parity states tt = 



• oscillator shells with N quanta are characterized by the constraint N = 2n + K. 

Thus for a given hyperangular and rotational configuration the quantumnumber n ladders 
the oscillator shells of increasing oscillator energy. 

2.1 Asymptotic solutions in coordinate representation 

The MJM method for solving the Schrodinger equation for quantum scattering systems is 
based on a matrix representation of the Schrodinger equation in terms of a square integrable 
basis. In a nuclear context the Harmonic Oscillator basis is appropriate. The boundary 
conditions are ultimately formulated in terms of the asymptotic behavior of the expansion 
coefficients of the wave function. 

In the case of three-cluster calculations, one needs to determine the asymptotic behavior 
of the wave function. Consider an expansion of the relative wave function 



with u = {n, K, {lil2)LM} and {^'y} a complete basis of six-dimensional oscillator states. It 
covers all possible types of relative motion between the three clusters. 

To obtain the asymptotic behavior of the three-cluster system the folding approximation 
is used. The assumption that antisymmetrization effects between clusters are absent in the 
asymptotic region is a natural one. It simplifies the many-body dynamics effectively to 
that of a three-particle system, as the cluster descriptions are frozen. The relative motion 
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problem of the three clusters in the absence of a potential (i.e. considering the JM reference 
Hamiltonian only) can then be explicitly solved in the HH method (see for instance [Hll^llUj). 
It involves the transformation of the Jacobi coordinates qi and q2 to the hyperradius p and 
a set of hyperangles Vl. The hyperspherical coordinates p and f2 describe the geometry of the 
three-particle system in the same way that the spherical coordinates describe the geometry 
of two-particle systems. The inter-cluster wave function in coordinate representation is 
expanded in HH's H^j^ (f2) where z/q has been chosen as a shorthand for {lil2)LM, and which 
are the generalization of the spherical harmonics Ylm {&■, 

In the absence of the Coulomb interaction this leads to a set of equations for the hyperra- 
dial asymptotic solutions, with the kinetic energy operator as the JM reference Hamiltonian 
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The solutions can be obtained analytically and are represented by a pair of Hankel 
functions for the ingoing and outgoing solutions: 
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where 
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One notices that these asymptotic solutions are independent of all quantum numbers uq, and 
are determined by the value of hypermomentum K only. In particular, different i^-channels 
are uncoupled. 

When charged clusters are considered the asymptotic reference Hamiltonian should con- 
sist of the kinetic energy and the Coulomb interaction: 
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The matrix ||/C|| is diagonal with matrix elements K{K + A), and ||Ze//||, the "effective 
charge", is off-diagonal in K and [hh)- The solution matrix ||7^(p)|| reflects the coupling 
of i^-channels. A standard approximation for solving these equations is to decouple the K- 
channels, by assuming that the off-diagonal matrix-elements of ||-Z^e//|| are sufficiently small: 



2m 



d^ 
dp^ 



5 d KiK + A)' 
p dp p2 



+ ^~e}Rk,.Ap) = o 



(12) 



The constants Z^ff depends on K and uq and all parameters of the many-body system under 
consideration. We will restrict ourselves to this decoupling approximation, but it is to be 
understood that its validity has to be checked for any specific three-cluster system. 
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Table 1: Two- and three-cluster asymptotic solution parameters 
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The asymptotic solutions then become 

RkL (p) = -^W±^v,K+2 {±2ikp) /pi (13) 
where W is the Whittaker function and rj is the well-known Sommerfeld parameter 

As ?7 is a function of li and I2 through the parameter ^e//, the asymptotic solutions 
will now be dependent on K and z/q. 

A two-cluster description leads to the well-known one-dimensional radial equation, with 
asymptotic solutions 

(p) = ^^±..,L+i (±2^^P) /P 
where p now is the inter-cluster radial distance and r] is the Sommerfeld parameter with 

One can easily combine the two- and three cluster asymptotics in a single representation 

as 

(p) = -^W^ir,,x {±2tkp) /p V (15) 

where the parameters C, A, a and r], differ for the two- and three-cluster channels, and are 
summarized in Table El 



2.2 Asymptotic solutions in oscillator representation 

The JM relies on an expansion in terms of oscillator functions, and the asymptotic behavior 
of the corresponding expansion coefficients c,y. We will restrict ourselves to the three-cluster 
situation; the two-cluster results can be readily obtained from and the results of the 
section 12. 1[ 

It was conjectured (see for instance |2j) that for very large values of the oscillator quantum 
number n the expansion coefficients for physically relevant wave-functions behave like 

Cn = {n\^) ~ V^pli^ibpn) (16) 
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with b the oscillator parameter of the basis, p„ = y/4n + 2K + 6 the classical turning point, 
and ip the hyperradial wave function. 

In the case of neutral clusters this leads after substitution of the hyperradial asymptotic 
solutions to the following expansion coefficients cit^ 

This result can be obtained in an alternative way ^21 by representing the Schrodinger 
equation, with the kinetic energy operator T as the reference Hamiltonian to describe the 
asymptotic situation, in a (hyperradial) oscillator representation 

oo 

J2 {n, {K, Uo)\f-E\m, {K, z/q)) c^'"' = (18) 

This matrix equation is of a three-diagonal form because of the properties of T and the 
oscillator basis. Solving for the expansion coefficients c^'*^" leads to a three-term recurrence 
relation 

rpK,iyo K,iyo , (rpK,UQ _ rp\ K,Uo , rpK,Uo K,Uo _ n /iqN 

n,n-l^n-l ' \'^n,n ^ j ' n,n+l^n+l ~ \^^) 

where 

T„^„^« = (n, (K, uo) |t| m, (K, z/q)) (20) 

The asymptotic solutions (i.e. for high n) of this recurrence relation are then precisely given 
by (Hlj). 

When the Coulomb interaction is present we again apply ()16|) to obtain 

J±)K _ / W^rl,^^ i^ikbpn) / \ 

In this case the oscillator representation of the Schrodinger equation is no longer of a tridi- 
agonal form, and cannot be solved analytically for the asymptotic solutions to corroborate 
this result. 

It should be noted that the above elaborations are valid for relatively small values of 
momentum k to which correspond sufficiently large values of the discrete hyperradius p„, 
which defines the asymptotic region. So for any value of k one can determine a minimum 
value for n to be safely in the asymptotic region. 

As we consider an asymptotic decoupling in the (i^, uq) quantumnumbers, one will deal 
with asymptotic channels characterized by the {K, z/q) values. So only in the internal (or 
interaction) region will states with different K and z/q be coupled by the short-range nuclear 
potential and the Coulomb potential. The three-cluster system can therefore be described 
by a coupled-channels approach, where the individual channels are characterized by a single 
/T- value, and we will henceforth refer to these channels as "i^-channels" . 
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2.3 Multi-channel JM equations 

In the current many-channel description of the JM formulation for three-cluster systems, the 
channels will be characterized by the a specific value of the set of quantumnumbers K, z/q, 
whereas the relative motion of clusters within the channel is connected to the oscillator index 
n. We will use K henceforth as an aggregate index for individual channels, and assume it 
represents all K, Uq quantum numbers. 

The Schrodinger equation can be cast in a matrix equation of the form 



E 

K',m 



H -E 



m, K' 







(22) 



As we will consider an ^-matrix formulation of the scattering problem, the expansion 
coefficients are rewritten as 

= + Sk.kcI-^'' - Sk.kc^^^^ (23) 

where, for each channel K, the Cn''^ are the so-called residual coefficients, the ci^'*'^ are 
the incoming and outgoing asymptotic coefficients. The matrix element Sk^k describes the 
coupling between the current channel K and the entrance channel Ki. 

As shown in [T31 [13 CHI , the c"it^^ satisfy the following system of equations 



oo 

E 

m=0 



n,K 



Hn-E 



m,K) c, 



Po "™,0 



(24) 



Hq being the asymptotic reference Hamiltonian, consisting of the kinetic energy operator for 
uncharged clusters, and the kinetic energy operator plus Coulomb interaction for charged 
clusters. The right-hand side for the equation for Cm^^ features /^q"'*'^ which is a regular- 
ization factor to account for the irregular behavior of the Cg ''^. This factor allows one to 
solve ()24|1 for all values of n. The value of I3q can be obtained for both reference Hamil- 
tonians (i.e. with or without Coulomb). We introduced (3^^^ = to keep the equations in 
a homogeneous form. 

The set of equations (j^^ for the asymptotic coefficients can then be solved numerically 
to different degrees of approximation depending on the requested precision. The c^it^^ have 
the desired asymptotic behavior (cfr eqs [T7| and EH) . 

Substitution of in the equations then leads to the following system of dynamical 
equations for the many-channel system: 



K'm 



n,K 



H-E 



m, K 



E 

K' 



>K,,K 



rj-)'^^' (25) 



where the coefficients Vn"^^^^ , defined in |2j, are given by 



KK' 



E 

m=0 



V 



(26) 
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This system of equations should be solved for both the residual coefficients Cn^^ and the 
S'-matrix elements Sk'k- 

To obtain an appropriate approximation to the exact solution of ()25|) . we consider an 
internal region corresponding to n < and an asymptotic region with n > N. The choice 

of N is such that one can expect the residual expansion coefficients < Cn'^^ 



to be sufficiently 

small in the asymptotic region. Under these assumptions (^3)) reduces to the following set 
of + 1 equations (n = O..A^): 



K',m<N 



H-E 



m, K') cT' - Y: Sk^k,V^^^^^' = -4-)^ 5.,o5k.k - Vt'^^^ (27) 



K' 



The total number of equations for an entrance channel Ki amounts to N^h * {N + 1), 
and solving the set of equations by traditional numerical linear algebra leads to N^h * 

N residual coefficients jclf'*''^; K = i^'min-.i^jnax; n = O..A^ — l| and Nch S'-matrix elements 

{SxiK', K = _ft'inin..i^jnax}- The sct of cquatious has to be solved for all N^h entrance channels 
labelled by Ki. 



2.4 Matrix elements and the Generating Function Method 

We only present the general principles for calculating matrix elements in a three-cluster basis. 
We will do so by considering the Generating Function Method. The two main quantities of 
interest to solve the JM equations are: the overlap matrix, and the Hamiltonian matrix. The 
former is of importance because of the proper normalization of the basis states. The latter 
is decomposed into the kinetic energy operator, the matrix elements of which are obtained 
mainly by group-theoretical considerations, the potential energy operator, which in our case 
will be chosen to be a local, semi-realistic, two-body interaction based on a superposition of 
Gaussians, and the Coulomb interaction. 

The basic principle of generating functions is well-known from mathematical physics. A 
generating function or generator state depends on a parameter, referred to as the generating 
coordinate, in such a way that an expansion with respect to that parameter yields basis states 
as expansion terms. Let us explain the principle using a familiar example: the single-particle 
translated Gaussian wave functions, which are appropriate for two-cluster descriptions 

0(r|R) =exp|-^r2 + V2 R-r- ^R^l (28) 

with the translation parameter R acting as generator coordinate. The choice of parametriza- 
tion of the generator coordinate influences the quantum numbers of the individual basis states 
that are generated. In a Cartesian parametrization R = [R^, Ry, R^) one generates the fa- 
miliar Cartesian (f)n^{Rx)(l)ny{Ry)(pnARz) oscillator states. With a radial parametrization 
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R = KR, (where the inverted hat stands for a unit vector) the expansion yields 



frlRl 



J2 ^"'^ 



r) 



(29) 



nl,m 



An underlying mathematical connection exists between such expansions, group representa- 
tion theory and coherent state analysis jTTl |^ . We exploit the generating function principle 
to facilitate the computation of matrix elements. The matrix element of any operator be- 
tween generating states is a function of the generating coordinates on the left and right 



X(R, R') 



(r|R) 



X 



(rlR': 



(30) 



Expansion of this function will yield the matrix elements between the basis states. They can 
be identified in the expansion by the appropriate dependence on the generator coordinates 



X 



■'n'l'm' 



(r)) (31) 



nlm n'l'm' 



For the three-cluster basis we consider the following generating function for inter-cluster 
basis functions (in what follows we shall use small q for the Jacobi vectors and capital Q for 
the corresponding generating coordinates) 



^ (qi, q2 IQi, Q2) = exp |-i (q? + q^) + ^2 (Qi ■ qi + Q2 ■ qs) " ^ (Q? + Q2) } 



(32) 



The choice of parametrization is linked to the basis states one intends to generate. Associated 
with our choice of basis (Zernike-Brinkman pi]), we introduce hyperspherical coordinates. 
The hyperradius and hyperangles, both for spatial coordinates and for generating parameters, 
are defined by: 



P = V + ?i = P cos e, q2 = p sin 9; 



R= ^Ql + Ql Q^ = RcosQ, Q2 = RsinG, (33) 
Using these, one expands the generating function (jH^ in HH functions: 

^(qi,q2|Qi,Q2) = 5^^.(p,^,qi,q2) s: (i?,e,Qi,Q2) (34) 

V 

where the full set of quantum numbers v (introduced previously) is involved in the summa- 
tion. The oscillator basis functions are 

vl/,(p,^,qi,q2) =Ar„,^ p^exp{V/2} ^f+'(/) i/jj^'^^^''' (0, qi, q2) (35) 



11 



and the generator coordinate functions are 



{R, e, Qi, Q2) = Mn,K H] 



r(hl2)LM 
K 



(0,Qi,Q2) 



(36) 



Here H denotes the HH function 



H 



r(hl2)LM 

K 



(e,Qi,Q2) 

^l'''^ (6) 



j^ihi2)LM^ii^i2) |Y,,(QO X Y,,(Q2)} 
(cos 9)'^ (sin 6)'^ PTJ:-\f' (cos 29) 



LM 



(37) 



From (jHBj) . one easily deduces the procedure for selecting basis functions with fixed quantum 
numbers v = {n, K, (/1/2) LM}. One has to differentiate the generating function {K + 2n)- 
times with respect to R and then to set i? = 0. After that one has to integrate over G with 
the weight to project onto the hypermomentum K; one has to integrate over unit 



The calculation of matrix elements with the generating function method is thus a two-step 
process. The first step is the calculation of the generating function for the operator involved. 
Usually this is accomplished with analytical techniques. The second step is the expansion 
of the generating function w.r.t. the generator coordinates. Either explicit differentiation or 
recurrence relations can be used to derive expressions for individual matrix elements In 
any case, the work involved is straightforward but extremely tedious; both approaches are 
best implemented using algebraic manipulation software such as Mathematica or Maple. 

A further explicit presentation of the calculation of matrix elements for overlap and 
Hamiltonian is beyond the scope of the current paper because of the highly technical details 
involved. We refer to |4j for explicit details for the three-cluster case. 



The numerical solution of the JM equations critically depends on an appropriate choice of 
N, separating the internal from the external region. The determining factor in this is the 
behavior of the potential energy matrix elements which can be slowly decreasing as a function 
of n. If this is the sufficiently large value of N has to be chosen. 

In the case of three-cluster systems it is known from literature f2I3 that the potential 
asymptotically behaves as in the hyperradius, with a corresponding effect on the matrix 
elements. It will be shown later that the asymptotic form of the effective potential in our 
calculation indeed follows this behavior. Potentials with an asymptotic tail dramatically 
affect the phase shift in the low energy region 122] and special care should be taken to 
get convergent results. This usually requires more than merely choosing a sufficiently large 
value of N. 




2.5 The Modified JM approach 
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Another factor that aversely affects convergence is the choice of a single basis (i.e. a 
fixed oscillator length b) for both the cluster and relative motion basis functions. This 
choice has been made to make the calculation of matrix elements manageable, and within 
acceptable execution time. The value chosen for b is usually fixed so that the frozen clusters 
have appropriate physical properties (e.g. ground-state energy). It is indeed important that 
the individual fragments are well described, as their properties cannot be influenced by the 
solution of the dynamical equations which only involves their relative behavior. This choice 
of basis is often far from optimal for a convergent calculation. 

To account for these convergence problems we have considered the MJM approach, de- 
veloped in In this approach we consider a semi-classical approximation for the matrix 
elements of the potential that is incorporated in the three-term recursion relation. We so 
introduce an intermediate region between the pure interaction and asymptotic regions. In 
this intermediate region a modified asymptotic solution holds, obtained from the modified 
recursion relation, and incorporating potential effects. In particular a long asymptotic tail 
of the potential can be taken into account in a satisfactory way. The matching point for 
the boundary condition, originally positioned at the border between interaction and asymp- 
totic regions, is brought down to the border between the near interaction and intermediate 
regions. This dramatically reduces the dimensions of the remaining set of matrix equations, 
as well as improves overall convergence by taking into account the asymptotic tail of the 
potential energy. 

3 Resonances in the MJM ^He and ^Be in the three- 
cluster approach 

In this section we consider the continuum obtained in a three-cluster MJM calculation for the 
6-particle nuclei ^He and ^Be, determined by the three-cluster configurations a + n + n and 
a + p + p. Our objective is to highlight the characteristics of MJM three-cluster calculations, 
and to produce accurate results for the astrophysically relevant resonances in ^He and ^Be. 

As the MJM for three-cluster scattering leads to a set of coupled-channel Schrodinger 
equations (see section ...), we will transform the 5- matrix to diagonal form. This is usually 
referred to as the eigenchannel representation of the ^-matrix. We will derive the position 
and width of the resonances eigenphase shifts 5 as a function of energy. To be precise, we 
use the conditions 

^,\e=e. = 0, r = 2 

We particularly focus on the properties of the low-lying 0"''-and 2"'"-resonances in ^Be and 
the 2"'"-resonance state in ^He. 

Several models and methods have already been applied to investigate the resonance states 
of ^He and ^Be. In the series of papers |1()| ESI 121], the nuclei were considered as three- 
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particle systems, neglecting antisymmetrization. These authors used an effective interaction 
between the alpha-particle and a nucleon, and a nucleon-nucleon (NN) realistic potential 
between the two nucleons. The HH's Method was used to describe the bound as well as the 
three-particle continuum states. In |2S] the Complex Scaling Method was used to calculate 
the characteristic parameters of the resonance states in ^He and ^Be. This was done within 
a three-cluster model, in which full antisymmetrization was taken into account, and the 
interactions between clusters was obtained from the Minnesota A^A^-potential. In [21] it was 
demonstrated that the method of continuation in the coupling constant (CCCM), used for 
the a + N + N configurations in ^He and ^Li, reproduces results, which are very close to 
those obtained by the CSM. We will compare MJM results to each of these methods. The 
MJM formulation leads to a good convergence of the computations, even for the (important) 
Coulomb contribution. An important renormalization effect on the irregular asymptotic 
solution still remains a restrictive factor in the reduction of the number of basis states in 
some cases. 

3.1 Details of the calculation 

The Volkov potential j2Z| has been used to describe A^A^ interaction in all calculations of 
this section. According to [HI IZ| this provides an acceptable description for ^He within the 
three-cluster model. The effective Volkov potential consists of central forces only, so that 
total angular momentum L, total spin S and total isospin T are good quantum numbers. We 
will therefore consider only three-cluster configurations a + A^ + A^ with S = and T = 1. 
This is known to be the most prominent spin-isospin state for the resonances of interest in 
the nuclei considered in this work. The Coulomb interaction has also been included, as it is 
to a great extent responsible for reproducing the O"*" resonance state in ^Be. 

The oscillator radius b, associated with the basis functions for the nuclear state, is the 
only free parameter for the MJM calculation. It's value was chosen to optimize the ground- 
state energy of the a-particle, and equals b = 1.37 fm. 

As explained before, we have a choice of Jacobi coordinate systems and this can be 
exploited to simplify the calculations. For a cluster configuration a -|- A^ -|- A^ one can 
consider the two Jacobi configurations displayed in figure ^ The first (referred to as the 
"4-1-2" -configuration) is the most appropriate for the current calculation. Selection rules 
significantly reduce the number of basis functions. Indeed, quantum numbers 5* = 0, T = 1 
coincide for the full six-nucleon system and the two-nucleon subsystem (A^ + A^). Hence the 
relative motion wave function must be an even function in the coordinate qi of this Jacobi 
system. This means that only even angular momenta /i have to be considered. Moreover, for 
positive parity states only even values of I2, and for negative parity states only odd values 
of I2 should be taken into account. With the second Jacobi configuration from figure ^ (the 
"5+1" one), these constraints are hard to meet and the full oscillator basis would have to 
be used. 

Table El shows the number of HH's {Nh) of given hypermomentum K for L'" = and 
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Table 2: Number (Nh) and accumulated number (Nc) of Hypersplierical Harmonics up to 
K = 10 



= 0+ 


K 





2 


4 


6 


8 


10 






1 


1 


2 


2 


3 


3 






1 


2 


4 


6 


9 


12 


L"" = 2+ 


K 





2 


4 


6 


8 


10 




Nh 




2 


3 


5 


6 


8 








2 


5 


10 


16 


24 



_ 2+ in the "4+2" Jacobi configuration. This table also indicates the total number of 
HH's so far {Nc) with K = K^i^, -^min + 2, . . . , -ft'max, i-e. it shows the number of channels 
for a given i^'max = K. 



3.1.1 Overlap matrix elements 

In the interaction region we apply full antisymmetrization, i.e. also inter-cluster anti- 
symmetrization. Here we look at its effect on positioning the boundary between the internal 
(interaction) region and the asymptotic region. No particle exchanges should occur in the 
asymptotic region. In the MJM the antisymmetrization effects can make themselves felt 
through the overlap and potential matrix elements. We use the notations of [4J for the over- 
lap matrix elements, and in particular the shorthand notation uq is used for the set {lil2)LM 
of quantum numbers: 

(n,(K,z/o)|l|n,(ir',z/o)) (39) 

Non-zero matrix elements ()39j) can be obtained from states within the same many-particle os- 
cillator shell only. As the oscillator shells in the Hyperspherical description are characterized 
by = 2n + i^' , the selection rule becomes 2n + K = 2n' + K' . 

In figure |21 overlap matrix elements diagonal in n for L = and hypermomenta K = 
and K = 2 are shown for the six-nucleon three-cluster system. One notices from this figure 
that the Pauli principle involves oscillator states of at least the 25 lowest shells, for both 
diagonal and off-diagonal matrix elements in K. The antisymmetrization effects are visible 
in the deviation from unity for the diagonal matrix elements, and the deviation from zero of 
the off-diagonal (in K) elements. These effects decrease monotonically with higher n. 

In figure El we compare matrix elements for L = diagonal in n and K for some of 
the i^'-values, with a a multiplicity quantum number for states with identical K. Only 
those matrix elements where the Pauli principle effect is most prominent have been shown. 
Some states with K = A and K = are affected more strongly by antisymmetrization than 
others. To understand this we note that the Hyperspherical angles (corresponding to the 
hyperangular quantum numbers hjh, LM) define the most probable triangular shape and 
orientation in space of the three-cluster system. The HH's with = 4, a = 2 (characterized 
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Figure 2: Matrix elements of the antisymmetrization operator for the O'^-state in ^He and 
^Be 
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Figure 3: Matrix elements of the antisymmetrization operator for the 0^-state in ^He and 
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Figure 4: Matrix elements of the antisymmetrization operator for the 0"'"-state of ^He and 
^Be in the three-cluster a + N + N {K = 0) and the two-cluster configuration a + 2N 



by h = h = 2) and K = 8,a = 3 (characterized by h = h = 4) seem to describe a triangular 
shape where one of the nucleons is very close to the a-particle. 

For larger i^- values the probability to find all clusters close to one another within a 
hypersphere of fixed radius p decreases, and one can expect that HH's with large values of 
K will play a diminishing role in the calculations. 

It interesting to compare overlap matrix elements for the three-cluster configuration 
a -|- -|- with those of the two-cluster configurations a + 2N (such as a -|- c? in ^Li or 
a + 2n'm ^He). This comparison is shown in figure |2 and it indicates that the Pauli principle 
has a much larger "range" in the three-cluster than in the two-cluster configuration. 

3.1.2 Potential matrix elements 

To investigate the Pauli effect on the potential energy of the three-cluster configuration 
a + N + N we compare the potential matrix elements with hypermomentum K = with full 
antisymmetrization against those in the folding approximation, where antisymmetrization 
between clusters is neglected. Figure |S] shows the diagonal potential matrix elements diagonal 
in n for K = and figure IHl shows the matrix elements along a fixed row n = 50 for K = 0. 
One notes that the folding model results are very close to the fully antisymmetrized ones, 
especially for larger n. 

In figure m we also display the potential matrix elements between states of the two lowest 
values of hypermomentum K = and K = 2. One notices that the K = 2 contribution 
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Figure 5: Diagonal potential matrix elements for the 0'*"-state of ^Be and ^He {K = 0) with 
full antisymmetrization, and in the folding model 
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Figure 7: Diagonal potential matrix elements for the 0"''-state of ^Be and ^He, with full 
antisymmetrization 



is the largest. The potential energy for = is relatively small, and so is the coupling 
between K = and K = 2 states. 

The main conclusion is that, in the asymptotic region, the "exact" potential energy 
can be substituted with the folding approximation. It leads to a considerable reduction in 
computational effort. We are led to the following setup for three-cluster calculations. In the 
internal region, consisting of states of the lower oscillator shells and with a large probability 
to find the clusters close to one another, the fully antisymmetrized potential energy is used. 
In the asymptotic region, where the average distance between clusters is large, we use the 
folding model potential. 

The folding model provides additional insight in the structure of the interaction matrix. 
It is well known that neither the a+n nor the n+n interaction can create a bound state in the 
corresponding subsystems of ^He. Only a full three-cluster configuration a + n + n contains 
the necessary conditions to create a bound state. In the folding model the total potential 
energy is a sum of contributions from the three interacting pairs: the a-particle and the first 
neutron, the a-particle and the second neutron and the neutron-neutron pair. For K = 
the first two contributions are identical, and we only have to consider the a + n and the 
n + n pairs. Figure |H1 shows the diagonal matrix element contribution of both components, 
and one notices that a + n represents the main contribution. 
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Figure 8: Folding model contributions in the diagonal potential matrix elements of the energy 
operator for the 0"^-state of ^Be and ^He {K = 0) 



Table 3: Effective charge matrix for the 0"*'-state in ^Be 



K- h,l2 


0;0,0 


2; 0,0 


4; 0,0 


4; 2, 2 


0;0,0 


7.274 


0.006 


-0.129 


1.414 


2; 0,0 


0.006 


7.146 


-0.436 


0.314 


4; 0,0 


-0.129 


-0.436 


7.428 


-0.877 


4; 2, 2 


1.414 


0.314 


-0.877 


9.098 



3.1.3 The effective charge. 

When Coulomb forces are taken into account, one needs to determine the effective charge 
in order to properly solve the MJM equations. The effective charge unambiguously defines 
the effective Coulomb interaction in each channel as well as the coupling between different 
channels. Using the approach suggested in section 12.11 the effective charges for the O"*"- and 
2+-states of ^Be were calculated. Part of the corresponding matrices of ||-^^H are displayed 
in tables El and m respectively. One notices that the diagonal matrix elements are much larger 
than the off-diagonal ones. This justifies the approximation ^ of disregarding the coupling 
of the channels in the asymptotic region. 

It is interesting to compare the three-cluster effective charge with the effective charge in 
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Table 4: Effective charge matrix for the 2"'"-state in ^Be 





2; 2,0 


2; 0,2 


4; 2,0 


4; 0,2 


4; 2, 2 


2; 2,0 


7.253 


0.400 


-0.224 


0.546 


-0.751 


2; 0,2 


0.400 


7.244 


-0.309 


-0.004 


-0.601 


4; 2,0 


-0.224 


-0.309 


6.942 


-0.186 


0.431 


4; 0,2 


0.546 


-0.004 


-0.186 


7.694 


-0.671 


4; 2, 2 


-0.751 


-0.601 


0.431 


-0.671 


7.345 



the two-cluster configuration. For the latter we can write 

Z = ZxZ^e^ 



2 / 



where Ai and Zi {A2 and Z2) are the respective mass and charge of both clusters. For the 
configuration a -|- 2p in ^Be, we then obtain an effective charge 

Z = 4=e^ - 6.65 
V3 

which is independent on the angular momentum of the system. One notice that the two- 
cluster effective charge is close to the diagonal matrix elements of ||-^^'||- We can assume 
that, if in one of the three-cluster channels the effective charge is very close to the two- 
cluster one, it could indicate that the two protons move as an aggregate in the asymptotic 
region. For the 2+-state we observe at least one channel with this property, carrying the 
labels X = 4, Zi = 2, Z2 = 0. 



3.2 Results 

3.2.1 Definition of the model space 

The model space for the current calculations is primarily determined by the total number of 
HH's in the internal and external region, and the number of oscillator states. Different sets 
of HH's can be used in the internal and asymptotic regions. An extensive set of HH's in the 
internal region will provide a well correlated description of the three-cluster system due to the 
coupling between states with different hypermomentum. The HH's in the asymptotic region, 
which are exactly (without Coulomb) or nearly exactly (Coulomb included) decoupled, are 
responsible for the richness in decay possibilities. 

In the current paper we restricted both internal and asymptotic HH sets maximal hy- 
permomentum value XmL = -f^max = 10. By extending the respective subspaces up to these 
maximal values, we obtain a fair indication of convergence. 

We fixed the matching point between internal and asymptotic region at = 50. It is 
based on the observations of the previous sections, i.e. (1) it is sufficiently large to allow the 
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E (MeV) 

Figure 9: Eigenphase shifts for the 0"'"-state of ^Be for different Ki 



Table 5: 0^-resonance parameters for ^Be, varying Kmlx {Kmax = 8) 



-limax 





2 


4 


6 


8 


E, MeV 


1.434 


1.314 


1.304 


1.298 


1.292 


r, MeV 


0.075 


0.082 


0.084 


0.085 


0.087 



Pauh principle its full impact, and (2) it is large enough for the semi-classical approximations 
in the MJM to be valid. 

The resonance parameters are obtained from the eigenphase shifts obtained from the 
eigenchannel representation (diagonalized form) of the S'-matrix. 

3.2.2 Convergence study 

We first consider the infiuence of the number of asymptotic channels. We do this by cal- 
culating the position and width of the 0~*"-state in ^Be using successively larger asymptotic 
subspaces. The first calculation includes all HH's up to iCmax = 8 and extends the asymp- 
totic subspace from K^lx = to Kmlx = 8. The second calculation includes all HH's up 
to K^nLx = 10 and extends the asymptotic subspace from K^lx = to Kmlx = 10. The 
corresponding results are shown in tables El and El One learns from these results that a 
sufficient rate of convergence has been obtained. Figure displays the first eigenphase shift 
as a function of energy for i^mL = 10 and a choice of Kmlx values from which results in the 
tables are derived. 
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Table 6: 0"'"-resonance parameters for ^Be, varying Kmlx {Kmax = 10) 



-limax 





2 


4 


6 


8 


10 


E, MeV 


1.324 


1.204 


1.192 


1.184 


1.176 


1.172 


r, MeV 


0.068 


0.069 


0.071 


0.071 


0.073 


0.072 



Table 7: 0+-resonance parameters for ^Be, varying Kmax {Kmax = 0) 



-''■max 





2 


4 


6 


8 


10 


E, MeV 




2.408 


2.020 


1.688 


1.434 


1.324 


r, MeV 




0.147 


0.129 


0.097 


0.075 


0.068 



Whereas the inclusion of higher hypermomenta in the asymptotics shows a fast and 
monotonia convergence, it is also clear from these results that a sufficient number of HH's 
has to be used for a correct description of the correlations in the internal state. To support 
this conclusion, we performed a calculation, again for the 0^-state in ^Be, in which only 
one HH with value K^"^ = was used. In the internal region the number of HH's was 
varied from K^^^ = up to if = 10. These results appear in table [Hand corroborate the 
previous conclusion. In particular they indicate that the effective potential obtained with 
the K^'^'^ = HH only, is unable to produce a resonance. Only after including a K = 2 
HH does the resonance appear. Further inclusion of higher HH states then lead to an 
acceptable convergence in a monotonically decreasing fashion for both position and width of 
the resonance. 

3.2.3 Comparisons 

In figure El we display eigenphase shifts for L'" = 0^ in ^Be in the full calculations, i.e. 
with the maximal number of internal and asymptotic HH's. One notices that the first 
-resonance state of ^Be appears in the first eigenchannel, and that a second (broad) 
resonance at a higher energy is created in the second eigenchannel (see table |H1 for details on 
this resonance). The phase shifts in the higher eigenchannels show a smooth behavior as a 



Table 8: Second resonance state parameters, obtained by MJM, HHM and GSM 



Nucleus 




MJM 


HH[23J 








E, MeV r, MeV 


E, MeV r, MeV 


E, MeV r, MeV 


^He 


0+ 


2.1 4.3 


5.0 6.0 


3.9 9.4 


^He 




3.7 5.0 


3.3 1.2 


2.5 4.7 


^Be 


ot 


3.5 6.1 






^Be 


2t 


5.2 5.6 
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Table 9: Resonance state parameters for ^He and ^Be, obtained by MJM, HHM, CSM and 
CCCM 





^He; = 2+ 


^Be; = 0+ 


^Be; = 2+ 


Method 


E, MeV r, MeV 


E, MeV r, MeV 


E, MeV r, MeV 


MJM 


1.490 0.168 


1.172 0.072 


3.100 0.798 


HHM [23 


0.75 0.04 






CSM [25J 


0.74 0.06 


1.52 0.16 


2.81 0.87 


CCCM^ 


0.73 0.07 







Table 10: Comparison of resonance state parameter in ^He and ^Be between MJM and 
experiment. 





MJM 


Experiment [30 




E, MeV r, MeV 


E, MeV r, MeV 


^He; L"" = 2+ 


1.490 0.168 


0.822±0.025 0.133±0.020 


^Be; = 0+ 


1.172 0.072 


1.371 0.092±0.006 


^Be; = 2+ 


3.100 0.798 


3.04±0.05 1.16±0.06 



function of energy without a trace of resonances in the energy range that we consider. 

Table IHl compares the results of this work to those obtained in other calculations, in 
particular by CSM by HHM |22 and by CCCM In tableUHlour results are compared 
with the experimental data available from [20] • The agreement with the experimental energy 
and width of the resonant states is reasonable. The difference between the experimental and 
calculated energies of the 2"'"-resonance states in both ^He and ^Be are probably due to the 
lack on LS-forces in the present calculations. 

It has been pointed out in and [23], that the barrier created by the three-cluster 
configuration is sufficiently high and wide to accommodate two resonance states, the second 
one usually being very broad. At first sight, this could be ascribed to an artifact of the 
HH Model since large values of hypermomentum K create a substantial centrifugal barrier. 
However the CSM calculations fIE\ I29j. which do not use HH's, also reveal such resonances. 
Our calculations now also confirm the existence of a second very broad resonance. A com- 
parison of the resonance parameters with those of the HH-method [23] and of the CSM [2^] 
is shown in table [Hi The differences are most probably due to the different descriptions of 
the system, as well as to the difference in iVA^-forces. 
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E (MeY) 

Figure 10: The eigenphases for the 0"'"-state in ^Be 



4 The ^H{^H,2nfHe and ^He{^He,2p)^He reactions in 
the MJM cluster approach 

A second apphcation that shows the strength of the MJM approach is connected to the 
nuclear reaction ^He{^He, 2p)'^He which contributes for 89% to the pp-chain of nuclear 
synthesis, and thus is of particular interest to astrophysics. We will couple the two- and 
three-cluster descriptions to calculate the reaction properties, by considering a two-cluster 
entrance, and a three-cluster exit channel 

The theoretical analysis of the ^He{^He, 2pYHe reaction is usually linked to its mirror 
companion, ^H{^H, 2n)^if e. A comparison of both leads to a better understanding of the 
underlying dynamics and of the Coulomb effects of reactions with three-cluster exit channels. 

A first microscopic calculation for these reactions was presented in A two-cluster 
approach was used for both the entrance and exit channels. The nucleon-nucleon fragment 
cluster (denoted A^A^ for either pp or nn) carried a simple shell-model description, featuring 
a pseudo-bound state with positive energy. The experimental cross-section or S'-factor at 
relatively high energy (approximately 1 MeV) was reproduced by adjusting the Majorana 
exchange parameter of the effective A^A^-potential. The available experimental data at the 
(small) energy range relevant in astrophysical reactions were fairly well reproduced. In 
this model no resonance state appears that would sufficiently amplify the S'-factor in the 
appropriate energy range, and thus would constitute an explanation for the solar neutrino 
problem. 
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Figure 11: Ground state of ^He as a function of the number of oscillator shells Ngt in the 
MJM three-cluster model, compared to the results of jS2]- The energy is relative to the 
a + n + n threshold. 



The previous model was further improved or enhanced in OH EH] by using a more 
elaborate description of the A^A^-channel, or by simulating the exit channel description using 
both the {'^He + A^) + and the ^ife + (A^ + A^) two-cluster configurations. The relative 
motion of the two clusters was described by a discrete superposition of translated Gaussian 
functions. In each case the same S'-factor shape as a function of energy was obtained. 

In this section we study both the {^H, 2n)^ He and ^He {^He, 2p)^ He reactions and 
use a correct treatment of the corresponding three-cluster exit channels. 

The bound state energy of ^He in the three-cluster description can easily be obtained 
by diagonalizing the nuclear Hamiltonian. We have compared it to the results obtained 
for the three-cluster calculation with the Stochastic Variational Method (SVM) \^2,^. 
We have used the Minnesota potential without spin-orbit components, and an oscillator 
parameter b = 1.285 fm which minimizes the a-particle energy as in In figure ITTl we 
compare our bound-state energy of ^He as a function of number of oscillator shells Nsh to 
the SVM value. We used all HH's up to Kmax = 10. One notices convergence for Ngh > 25 
towards E = —0.8038 MeV (relative to the a + n + n threshold). This is to be compared 
to E = —1.016 MeV for the full SVM calculation; full convergence towards the SVM result 
would require additional K values, which is beyond the scope of the current calculation. 
These results encourage us to combine the two- and three-cluster MJM descriptions to obtain 
an advanced description of the fusion reactions ^H{^H, 2nYHe and ^Hei^He^ 2pYHe. 
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4.1 Model specifics 



We again rely on section |2l for the details of the microscopic model. The specific cluster 
configurations used to describe the six-nucleon systems ^He and ^Be are analogous to those 
of section El 

4.1.1 A combined cluster model 

The six-nucleon wave functions will be built up by using both the two- and three-cluster 
configurations, each one fully antisymmetrized: 

mL = A{^^N ^37V /L(qo)}+^{^a N ^L(qi,q2)} (40) 

The \E'a {N stands for either nucleon, a for ^ife, and 3A^ for '^H or ^He) represents cluster 
component wave functions, fi and ql refer to the wave functions of relative motion for the 
two-, and three-cluster system. The are Jacobi coordinates describing the configuration 
of relative position of the clusters. 

All of the reaction dynamics, and in particular the S'-matrix elements, is concentrated 
in the functions describing the relative motion, i.e. fi and ql because the internal cluster 
wave functions are "frozen". 

We need to make a choice for the Jacobi coordinates q and for a classification scheme of 
the wave-functions to be used in the expansion of fi and ql- 

For the two-cluster configuration [ST] (in '^H +^ H respectively ^He +^ He) we use the 
standard spherical coordinates qo = {q'o, qo} and take the quantum numbers n = {n, L, M} 
to classify the basis states. The n is the radial oscillator quantum number. As we have only 
central components in the NN interaction, the angular momentum L of relative motion will 
be an integral of motion for the system. 

For the three-cluster configurations (in ^He + p + p respectively ^He + n + n) we use the 
Hyperspherical coordinates (jSS|) . This choice is consistent with the set of quantum numbers 
z/ = {N,K,{lil2)LM} = {np,K,{lil2)LM}, in which N = 2np + K represents the total 
number of oscillator quanta, and Hp reflects the number of hyperradial excitations. The 
partial angular momenta li and I2 are associated with the choice of Jacobi vectors qi and 
q2. A coupled ii'-channel calculation with each channel characterized by the set of quantum 
numbers z/q = {K, {lil2)LM}, has to be performed. This type of basis is particularly suitable 
for the so-called Borromian nuclei, and nuclei with pronounced three-cluster features, when 
the three-cluster threshold represents the lowest energy decay channel. 

4.1.2 The boundary conditions 

The MJM boundary conditions are expressed in terms of the expansion coefficients of the 
wave functions of relative motion. They are directly connected to the boundary conditions 
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in coordinate representation. For the two-cluster configurations the asymptotic form of the 
expansion coefficients in = ^ Cn,L(pn,L can be approximated by 

C„,L ^ aA^/l (r„) (41) 



where the {4>n,L} are the oscillator basis functions, and r„ = by/ 4n + 2L + 3 is the classical 
turning point of the three-dimensional oscillator with energy En = hw {2n + L + 3/2). In 
three-cluster configurations the expansion is defined hj ql = J2 dnp,L4>np,L, where the {(pnp,L} 
now stand for the oscillator basis for three-cluster relative motion. The expansion coefficients 
behave asymptotically as 

rf„^,i ^ V2pl9L (Pn) (42) 



with Pn = h^jAup + 2K + 6. For clarity we have indicated in the preceding discussion only 
relevant indices to this section. 

We consider both incoming and outgoing waves for the two-cluster configurations 

fL (qo) ^ [4'^ (^ogo) - S{,},{,} (fcogo)] Ylm (qo) (43) 

where S^^-^^i^^} is a notation to characterize the elastic two-cluster scattering matrix element 
for the +^ H, resp. ^He +^ He channels, and Y^m (qo) is the spherical harmonic. 

Because we are only interested in reactions with a three-cluster exit-channel, the asymp- 
totic wave function can be written as 

9L (qi, q2) = 9L (p,^) ^ Yl [-^M^M (kp)] H''^ (fi) (44) 

where S^^}^^^^^ is the scattering matrix element describing the inelastic coupling between the 
two- and three-cluster channels, and (fi) is the HH. 
The total cross-section is given by the expression 

vr (2L + 1)(25 + 1) ,2 ,,,, 

^ = P E ^ T E l%}>o}| (45) 

° L,S UQ 

with S the total spin of the six-nucleon system. 

As discussed in section 12.11 the asymptotic solutions for incoming and outgoing waves 
can be written as 

{kp) = -^W±,,,x i±2tkp) /p^ (46) 

where W are the Whittaker functions, and r] the Sommerfeld parameter (for the parameters 
C, A, cr and rj for two- and three-cluster channels: see table |2I). 

By using the correspondences (jlT|l and we can now define the boundary conditions 
for the expansion coefficients 



Cn,L 



,(+) 



dnp,uo ^ Pn -%},{^o} '^K ikpn) (47) 



28 



or equivalently 



using the notations 

<,lo ^ Pn ^i?^ (^Pn) (49) 

The matching of internal and asymptotic regions is equivalent to the one in the tradi- 
tional Resonating Group Method (RGM). The correspondence between the matching point 
in coordinate space for RGM and in function space for MJM is easily made (see through 
the value of the classical oscillator turning point r„ = b\^An + 2L + 3 for 2-cluster systems 
and pn = b^y^Up + 2K + 6 for 3-cluster systems. An appropriate value for the matching 
point can be obtained by choosing sufficiently large values for the total number of oscillator 
quanta N = 2n + L = 2np + i^' in the internal region. 



4.1.3 Shape analysis 

The HH's can reveal information on the spatial distribution of clusters and of the reaction 
dynamics. 

These harmonics define a probability distribution in 5-dimensional coordinate (momen- 
tum) space for fixed values of hyperradius: 

dW^^ ( Q) = |n„ ( Q)f dQ, dW^^ ( Qk) = 1^0 ( ^k)f dQk (50) 

By analyzing the probability distribution, one can retrieve the most probable shape of three- 
cluster shape or "triangle" of clusters. A full analysis of a function of 5 variables is non-trivial 
and one usually restricts oneself to some specific variable(s). We integrate the probability 
distribution dW^^ ( Q) over the unit vectors qi,q2 (resp. ki,k2) 

dW^, (9) = j cos^ 9 sin^ 9d9 rfqirfqz 

dW,, {9k) = j |n„ ( CDs'" 9ksm^ 9kd9k rfkirfka (51) 
and introduce the (new) variable(s) 

£ = 4 = cos2^, 8 = % = cos'9, 
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Table 11: Number of Hyperspherical Harmonics for L = 0. 
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In coordinate space these can be interpreted as the squared distance between the pair of 
clusters associated with coordinate qi, or, in momentum space, the relative energy of that 
pair of clusters. We obtain 



dW,, {d) 



dO 



N 



ih,h) 
K 



cos'^ ^sin'^ ^pfe+i/2A+i/2) (cos 2^) 



COS 6' sin 6 



N 



K 



'1 



gy^/"^ p(/2+l/2,Zi+l/2) 



(2^-1) V^(l-^) 



(52) 



This function represents the probability distribution for relative distance between the two 
clusters, resp. for the energy of relative motion of the two clusters. The kinematical factor 
cos^ 9 sin^ 9 was included to make Wi,^ (S) proportional to the differential cross section in 
momentum space, provided the exit channel is described by the single HH Y,yg (Q). 

In figure El we display W^o (S) for some HH's involved in our calculations. These figures 
show that different HH's account for different shapes of the three-cluster systems. For 
instance, the HH with K = 10 and li = I2 = prefers the two clusters to move with very 
small or very large relative energy, or, in coordinate space, prefers them to be close to each 
other, or far apart. 



4.2 Results 

Again we use the VP as the NN interaction. The Majorana exchange parameter m was 
set to be 0.54 which is comparable to the one used in [33]. The oscillator radius was set to 
b = 1.37 fm (as in [HIIZ|) to optimize the ground state energy of the alpha-particle. 

The VP does not contain spin-orbital or tensor components so that total angular momen- 
tum L and total spin S are good quantum numbers. Moreover, due to the specific features 
of the potential, the binary channel is uncoupled from the three-cluster channel when the 
total spin 5* equals 1; this means that odd parity states L'" = 1", 2~, . . . will not contribute 
to the reactions. 

To describe the continuum of the three-cluster configurations we considered all HH's with 
K < Kmax = 10. In Table ITD we enumerate all contributing i^-channels for L = 0. For each 
two- and three-cluster channel we used the same number n = Hp = Ni^t of basis functions 
to describe the internal part of the wave function \E'i. Ni^t then also defines the matching 
point between the internal and asymptotic part of the wave function. We used Ni^t as a 
variational parameter and varied it between 20 and 75, which corresponds to a variation in 
coordinate space of the RGM matching radius approximately between 14 and 25 fm. This 
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Figure 13: 5'-factor of the reaction {^H, 2n) He. Experimental data are taken from 
|SH](Serov), pi (Govorov),[in] (Brown) and [IT] (Agnew) 



variation showed only small changes in the S-matrix elements, of the order of one percent or 
less, and do not influence any of the physical conclusions. We have then used Mint = 25 for 
the final calculations as a compromise between convergence and computational effort. We 
also checked the impact of Nint on the unitarity conditions of the S'-matrix, for instance the 
relation 

I'D 

We have established that from Nint = 15 on this unitarity requirement is satisfied with a 
precision of one percent or better. In our calculations, with Nint = 25, unitarity was never a 
problem. It should be noted that our results concerning the convergence for the three-cluster 
system with a restricted basis of oscillator functions agree with those of Papp et al [HZ|, where 
a different type of square-integrable functions was used for three-cluster Coulombic systems. 

In figure El we show the total S'-factor for the reaction '^H {^H, Inf He in the energy 
range < < 200 keV. One notices that the theoretical curve is very close to the experi- 
mental data. The total S'-factor for the reaction ^He {^He, 2p)^ He is displayed in figure IT^ 
It is also close to the available experimental data. The S'-factor for both reactions is seen to 
be a monotonic function of energy, and does not manifest any irregularities to be ascribed 
to a hidden resonance. Thus no indications are found towards explaining the solar neutrino 
problem. 
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Figure 14: 5'-factor of the reaction ^He {^He, 2p) He. Experimental data are from |l2] 
(Krauss), gSj (LUNA 99) and 44 (LUNA 98). 



The astrophysical ^-factor at small energy is usually written as 

S (E) = So + S'^E + S'^E^ (53) 

We have fitted the calculated S-factor to this formula in the energy range Q < E < 200 keV. 
For the reaction {^H, 2nf He we obtain the approximate expression: 

5 (^) = 206.51 - 0.53 E + 0.001 ^2 keV b (54) 

and for ^He {^He, 2pf' He we find: 

5(E) = 4.89 - 3.99 E + 2.3 10"^ ^2 MeV b (55) 

One notices significant differences in the S'-factor for the ^He and ^Be systems. The A^A^- 
interaction induces the same coupling between the clusters of entrance and exit channels 
for both ^He and ^Be. It is the Coulomb interaction that distinguishes both systems, and 
accounts for the pronounced differences in the cross-sections and S'-factors. 

We compare the calculated S'-factor to fits of experimental results for the reaction 
^He {^He, 2pf He: 



S{E) = 5.2-2.8 E + 1.2 E^ MeV b gS] 

S (E) = (5.40 ± 0.05) - (4.1 ± 0.5) E + (2.3 ± 0.5) E^ MeV bgHl 

S(E) = (5.32 ±0.08) - (3.7 ±0.6) (1.95 ±0.5)^2 MeV b gH (56) 
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Figure 15: Three-cluster channel contributions to the total S'-factor for the reactions 
{^H, 2n) He and ^He {^He, 2p) He in a full calculation with Kmax = 10. 




Figure 16: Three-cluster channel contributions to the total S'-factor of the reactions 
^H{^H, 2n)'^He in a full calculation with Kmax = 10, in the energy range < E < 1000 keV. 
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Figure 17: Convergence of the S'-factor of the reaction ^H{^H, 2nYHe for K^ax ranging 
from to 10. 



The constant and hnear terms of the fit display a good agreement. The difference in energy 
ranges between the calculated (0 < < 200 keV) and experimental (0 < < 1000 keV) 
fits make it difficult to attribute any significant interpretation to the discrepancy in the 
quadratic term. 

The HH's method now allows to study some details of the dynamics of the reactions con- 
sidered. In Figs. EJand^lwe show the different three-cluster i^-channel contributions {Wy^-^) 
to the total S'-factor of the reactions. In figure El these contributions (in % with respect to 
the total S-factor) are displayed for some fixed energy (1 keV), while figure UHl shows the de- 
pendency of Wy^ (in absolute value) on the energy of the entrance channel. One notices that 
three HH's dominate the full result, namely the {K = 0; li = I2 = 0}, {K = 2; li = I2 = 0} 
and {K = 4; li = I2 = 2}, and this is true in both reactions. The contribution of these states 
to the S-factor is more then 95 % . There also is a small difference between the reactions 
"^H {^H, 2n)^ He and ^He{^He, 2p)^i7e, which is completely due to the Coulomb interaction. 

The Figs. El and El yield an impression of the convergence of the results. We notice 
that the contribution of the HH's with K > 6 is small compared to the dominant ones. 
This is corroborated in figure El where we show the rate of convergence of the S'-factor in 
calculations with K^ax ranging from up to 10. Our full K^ax = 10 basis is seen to be 
sufficiently extensive to account for the proper rearrangement of two-cluster configurations 
into a three-cluster one, as the differences between results becomes increasingly smaller. 

To emphasize the importance for a correct three-cluster exit-channel description, we 
compare the present calculations to those in [3T], where only two-cluster configurations 
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Figure 18: Comparison of the S'-factor of the reaction ^H{^H, 2n)^He in a calculation with 
a three-cluster exit-channel and a pure two-cluster model. 

'^He + 2n resp. "^He + 2p were used to model the exit channels. In both calculations we used 
the same interaction and value for the oscillator radius. In figure ITHl we compare both results 
for ^H{^H, 2n)'^He. An analogous picture is obtained for the reaction ^He{^He, 2p)'^He. 

4.2.1 Cross sections 

Having calculated the S-matrix elements, we can now easily obtain the total and differential 
cross sections. In this section we will calculate and analyze one-fold differential cross sections, 
which define the probability for a selected pair of clusters to be detected with a fixed energy 
Ei2. To do so we shall consider a specific choice of Jacobi coordinates in which the first 
Jacobi vector qi is connected to the distance between these clusters, and the modulus of 
vector ki is the square root of relative energy £'12. With this definition of variables, the cross 
section is 



2 



da {E12) ~ — 




rfkidka S{^}{^^}Y^o (fife) sin^ 9k cos^ 9kd9k 



(57) 



After integration over the unit vectors and substitution of sin6'fc, cos9k, d9k with 
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Figure 19: Partial differential cross sections of the reactions ^H{^H, 2n)^He and 
^He{^He,2p)^He. 



37 



4 8 12 16 

^,2 (McV) 

Figure 20: Calculated and experimental differential cross section for the reaction 
^He{^He,2p)'^He. Experimental data is taken from ^45j . 
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one can easily obtains da {E12) /dE^. 

In figure^lwe display the partial differential cross sections of the reactions '^H {^H, 2rif He 
and ^He i^He, 2pf He for the energy E = IQ keV in the entrance channel. The solid lines 
correspond to the case of two neutrons (protons) with relative energy £^12, while the dashed 
lines represent the cross sections of the a-particle and one of the neutrons (protons) with 
relative energy £"12. 

We wish to emphasize the cross section in which two neutrons or two protons are simulta- 
neously detected. One notices a pronounced peak in the cross section around E12 ~ 0.5 MeV. 
This peak is even more pronounced for the reaction ^He {^He, 2p)^ He. It means that at such 
energy two neutrons or two protons could be detected simultaneously with large probability. 
We believe that this peak can explain the relative success of a two-cluster description for the 
exit channels at that energy. The pseudo-bound states of nn- or pp-subsystems used in this 
type of calculation then allows for a reasonable approximation of the astrophysical S'-factor. 

Special attention should be paid to the energy range 1-3 MeV in the ^He + n and ^He + p 
subsystems. This region includes 3/2~ and l/2~ resonance states of these subsystems with 
the Volkov potential. In figure ITIH (dashed lines) we see that it yields a small contribution to 
the cross sections of the reactions ^H {^H, 2nf He and ^He {^He, 2p)'^ He. This contradicts 
the conclusions of j^l] and where the 1/2^ state of the '^He + N subsystem played a 
dominant role. We suspect this dominance to be due to the interplay of two factors: the 
weak coupling between incoming and outgoing channels, and the spin-orbit interaction. 

In figureEOlwe compare our results for the total proton yield (reaction ^He {^He, 2p) He) 
to the experimental data from jlHl- The latter were obtained for incident energy E i^He) = 
0.19 MeV. One notices a qualitative agreement between the calculated and experimental 
data. 

The cross sections, displayed in Figs. ^1 and 1201 were obtained with the maximal number 
of HH's {K < 10). These figures should now be compared to the figure El which displays 
partial differential cross sections for a single ii'-channel. The cross sections, displayed in 
Figs. El and EDI differ considerably from those in Figs. El and comparable ones, even for 
those HH's which dominate the wave functions of the exit channel. An analysis of the cross 
section shows that the interference between the most dominant HH's strongly influences 
the cross-section behavior. To support this statement we display the proton cross sections 
obtained with hypermomenta K = 0, K = 2, K = 4 to those obtained with the full set of 
most important components K^ax < 4 in figure |^ One observes a huge bump around 10 
MeV which is entirely due the interference of the different HH components. We also included 
the full calculation [K^ax < 10) to indicate the rate of convergence for this cross- sect ion. 
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Figure 21: Partial cross sections for the reaction ^He(^He,2p)'^He obtained for individual 
K = 0,2 and 4 components, compared to the coupled calculation with K^ax = 4 and the 
full calculations with K^ax = 10. 

5 Conclusion 

In this paper we have presented a three-cluster description of light nuclei on the basis of 
the Modified J-Matrix method (MJM). Key steps in the MJM calculation of phase shifts 
and cross sections have been analyzed, in particular the issue of convergence. Results have 
been reported for ^He and ^Be. They compare favorably to available experimental data. 
We have also reported results for coupled two- and three-cluster MJM calculations for the 
^ife {'^He, He and {^H, 2nf He reactions with three-way disintegration. Again com- 
parison indicates good agreement with other calculations and with experimental data. 
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